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The flow past inline oscillating rectangular cylinders is studied numerically at a Reynolds number 
representative of two-dimensional flow. A symmetric mode, known as S-II, consisting of a pair of 
oppositely-signed vortices on each side, observed recently in experiments, is obtained computation- 
ally. A new symmetric mode, named here as S-III, is also found. At low oscillation amplitudes, 
the vortex shedding pattern transitions from antisymmetric to symmetric smoothly via a regime of 
intermediate phase. At higher amplitudes, this intermediate regime is chaotic. The finding of chaos 
extends and complements the recent work of Perdikaris et al. [1 . Moreover it shows that the chaos 
results from a competition between antisymmetric and symmetric shedding modes. Rectangular 
cylinders rather than square are seen to facilitate these observations. A global, and very reliable, 
measure is used to establish the existence of chaos. 



I. INTRODUCTION 



Vortex shedding from bluff bodies is an extensively 
studied problem. The preferred mode of vortex 
shedding, in the uniform flow past a fixed body, is 
antisymmetric. On the other hand a body forced 
to oscillate in quiescent fluid would be expected to 
shed a symmetric pattern of vortices. Thus, for a 
body oscillating inline in a uniform external flow, 
as the frequency of oscillation is increased with all 
other parameters remaining fixed, we may expect a 
transition from antisymmetric to symmetric shed- 
ding. Both kinds of shedding have been observed 
[9l [121 US] • Secondly, it has been seen experimen- 
tally that there is more than one kind of symmetric 
shedding [8 j. Our objective here is to improve our 
understanding of the frequency response of the sys- 
tem, in terms of the spatial arrangements of vor- 
tices and the transitions therein. Studies on such 
flows have found application in predicting the load- 
ing on offshore structures [5|. Also, this is a simple 
example of the flow due to an accelerating body. 
The prediction of flow patterns in its wake can be 
important in various contexts, such as in the track- 
ing of underwater bodies. 

Griffin & Ramberg [13] were among the first to 
study vortex shedding from an inline oscillating cir- 
cular cylinder in a freestream. They found that 
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the vortex shedding frequency locks on to the fre- 
quency of cylinder oscillation for 1.2 < f e /f Q < 2.5, 
where f e is the frequency of cylinder oscillation and 
f would have been the frequency of vortex shed- 
ding if the cylinder were held stationary. The sub- 
scripts e and o have been chosen to stand for 'exci- 
tation' and 'original' respectively. Both primary 
lock-on, where the shedding frequency f s = / e , 
and subharmonic lock-on, with f s = f e /2 were ob- 
served. Ongoren & Rockwell [12] carried out ex- 
periments with a circular cylinder oscillating at 
an angle a with a uniform freestream. Outside 
the lock-on regime, competition between symmet- 
ric and antisymmetric modes in the form of switch- 
ing of modes in a single experiment were observed. 
In some recent experiments Konstantinidis & Bala- 
bani [9 J too noted the symmetric mode mentioned 
above, where all vortices shed from the top wall 
were of one sign, while those shed from the bot- 
tom wall were of the opposite sign. This pattern 
is called the S-I mode of shedding. In another ex- 
perimental study on a circular cylinder, Xu et al 
[8] discovered a new mode of symmetric shedding, 
which they named S-II. Two vortices of opposite 
sense were shed from each side (top and bottom) 
during each cycle. This mode was observed for high 
frequencies and amplitudes. There was consider- 
able reverse flow during a part of the cycle, which 
aided in the formation of opposite signed vortices 
on a given side of the cylinder. Very few numerical 
studies have reported the symmetric S-I shedding, 
Zhou & Graham's [15] being one. To our knowl- 
edge, the S-II mode has not been found numerically 
before. 

Besides systematic shedding, we could have chaotic 
shedding. Chaos in flow around an inline oscillat- 
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ing circular cylinder (or equivalently, in oscillating 
flow past a fixed cylinder) was reported by Vittori 
& Blondeaux [2! and Perdikaris et al. [1 . The 
former study had no mean flow, and showed that 
the route to chaos is quasiperiodic, and the lat- 
ter study attributed chaos to mode competition. 
However no evidence of mode competition was pro- 
vided. One goal of the present study is to present 
direct evidence of competition between antisym- 
metric and symmetric shedding, and the resulting 
chaos. Ciliberto & Gollub [3j 0] showed that com- 
petition between different modes in parametrically 
forced surface waves can result in chaos. Their 
study also revealed the existence of chaotic 'win- 
dows'. They remarked that this could be a com- 
mon cause of chaos in systems in which different 
spatial structures can exist. The present flow is 
shown to be one example. 

In a circular or square cylinder, the behaviour of 
each row of shed vortices can be clouded by inter- 
action with the opposite row. For this reason, we 
study the flow past rectangular cylinders of vari- 
ous aspect ratios. To our knowledge, there have 
been no previous studies on vortex shedding from 
inline oscillating rectangular cylinders. We show 
the existence of symmetric S-I and S-II modes, and 
a new S-III mode, besides the Couder-Basdevant 
mode. We then discuss a physical mechanism for 
the enhancement of S-II shedding from rectangu- 
lar cylinders. We show that the transition from 
antisymmetric to symmetric shedding, as f e is in- 
creased, for lower oscillation amplitudes occurs via 
periodic flows of different phase. At higher oscilla- 
tion amplitudes we find windows of chaos between 
the regimes where shedding modes are antisymmet- 
ric and symmetric. 



time between two successive collisions. Thus the 
density p = and the momentum vector is 

pUk = TiifiCik- No-slip on walls is imposed by em- 
ploying the bounce-back scheme. 

It is straightforward to show that in the absence of 
rotation of the cylinder, the vorticity and continu- 
ity equations in a cylinder-fixed frame are the same 
as those in the lab-fixed frame. In the cylinder- 
fixed coordinate system we have an oscillating inlet 
flow [7 , u x = Uoo + AQ sin Fit, A is the amplitude 
of the displacement of the body, Q (= 2irf e ) is the 
frequency of oscillation, u x and u y are the stream- 
wise and transverse components of fluid velocity 
respectively. The other boundary conditions are 
du x /dx — duy/dx = at the outlet; u x = u y = 
on the cylinder; u x = + AQ sin Qt, u y = at 
the top and bottom surfaces. The characteristic 
length in this study is taken to be the height D of 
the body, and the characteristic velocity is Uqq, so 
the Reynolds number is defined as Re = UoqD/v, 
and the Strouhal number as St = f s D/Uoo, f s be- 
ing the shedding frequency of the vortices. 



Figure 1(a) shows a part of the domain used for the 



simulations. The size of the domain is 75D x 75D 
with nonuniform grids of up to a million grid 
points, whose arrangement is shown in figure [T(b)| 
The values for the domain size and number of grid 
points were arrived at after doing a grid indepen- 
dence study. The reason for the use of a large do- 
main is twofold: to minimize the effects of out- 
let boundary conditions on the wake and to ensure 
that the symmetric modes are insensitive to the lat- 
eral extent of the domain. All simulations in this 
study were carried out at Re = 200. 



II. PROBLEM FORMULATION 



The numerical procedure described in Dixit & 
Babu [16] was employed after making suitable 
changes. Flow solutions are obtained by solving 
the Lattice-Boltzmann equation 
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As in the usual nomenclature, /i(X, t) and 
f^ q (X,t) are the nonequilibrium and equilibrium 
distribution functions in i th direction, r is the 



The numerical approach was validated as follows. 
Strouhal numbers obtained from flow past a fixed 
square cylinder were found to be in good agree- 
ment for a range of Reynolds numbers from 50 to 
250 with those from Okajima's experiments [11 
and Ansumali et aZ.'s numerical simulations [T4] . 
With an oscillating square cylinder, the dominant 
frequencies in the spectrum for the wall normal ve- 
locity u y at a typical location obtained from our 
simulations was found to be in good agreement 
with the spectrum for the lift coefficient obtained 
by Minewitsch et al. [7]. Further details are avail- 
able in Srikanth et al. [10 . Next, as done in [7], 
the frequency ratio was fixed at f e /f Q — 1-6 and 
A/D was varied from 0.15 - 0.4 in steps of 0.05. In 
excellent agreement with [7J, the lock-on window 
between 0.15 < A/D < 0.4, where f s = 0.5/ e is 
reproduced. 

For Re = 200, A/D = 0.175 and f e /f Q = 2 in 
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FIG. 1: (a) Schematic of the domain used in the simula- 
tions. Time signals are stored at the monitor points A, 
B and C, located at (2/i,0), (h, D) and (2.25D, 0.5D) 
respectively (not to scale), (b) The grid used. 



(b) 

FIG. 2: Qualitative comparison of shedding modes for 
the square cylinder. Vorticity fields are shown, (a) 
f e /f Q = 2 & A/D = 0.175. This mode is similar to the 
one obtained by Couder & Basdevant [22]. (b) f e /fo = 
1.73 & A/D = 0.5. This is S-II mode of shedding, 
similar to that seen in the experiment of Xu et al. [8] . 



III. RESULTS 



A. Modes of vortex shedding 



the square cylinder case, the vorticity field qualita- 
tively matches with one of the experimental results 
of Couder & Basdevant [22]. This mode, shown in 



figure 2(a), consists of two rows: one with binary 
vortices and the other with single vortices. We are 
able to obtain the S-II mode o f vorte x shedding be- 
hind a square cylinder, figure |2(b)| for Re = 200, 
A/D = 0.5 and f e /f Q = 1.73. Xu et al [8 were 
the first to report the S-II mode of shedding ex- 
perimentally in the case of a circular cylinder. The 
amplitude ratio and frequency were the same in 
their experiment as the values used here, but their 
Reynolds number was 500. To the best of our 
knowledge, the S-II and Couder-Basdevant modes 
have not been seen computationally before. 



When the oscillation frequency f e is very small, the 
flow is not too different from that past a fixed cylin- 
der, except that the Reynolds number now is slowly 
varying. One expects, and finds, a slightly modi- 
fied Karman street behind the body in this case. 
The same is true when the amplitude of oscilla- 
tion A/D is small, since the oncoming flow merely 
sees a slightly modified body on an average. To ob- 
serve competition between symmetric and antisym- 
metric shedding, one needs an effective oscillation 
Reynolds number Re Q = (QA)D/v which is not 
negligible compared to that of the incoming flow. 
Although we carried out many simulations to en- 
sure that our results are general, we present only a 
few typical ones. The shedding pattern changes as 
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the excitation frequency f e is increased from 0.5/ o , 
half the natural shedding frequency of a stationary 
cylinder, to five times this value, some examples are 
shown in figure [3J A rectangular cylinder of aspect 
ratio 4 is used here and A/D is fixed at 0.1. At 
low f e the shedding is antisymmetric, and goes to 
symmetric shedding as f e increases. At / e // Q = 5 
we have the symmetric S-I mode, with all the top 
vortices being of one sign, and all the bottom of 
the opposite sign. At moderate f e /fo, the shed- 
ding is neither symmetric nor antisymmetric, but 
the upper and lower vortices are shed with a phase 
between and tt (or tt and 2tt). The flow how- 
ever is still periodic. In some cases, vortex merger 
on each side of the cylinder is promoted, and the 
pattern downstream becomes antisymmetric. At 
small / e , the shedding frequency f s is close to f Q . 
However, as f e /fo is increased beyond 2, f s de- 
creases before locking on to a subharmonic of / e , 
and then starts increasing proportionately with / e , 
such that f s /f e = 0.25 for 3 < f e /f Q < 4. Be- 
yond this range f s steadily decreases with further 
increase in f e . The lock-on is similar to those seen 
on circular cylinders [9]. Note that at these high 
frequencies shedding occurs on a given surface once 
every four complete oscillations, rather than once 
in every other oscillation. 

Next, choosing A/D = 0.175, we summarise in fig- 
ure [4] the patterns of vortex shedding observed on 
cylinders with aspect ratio D/h = 1,2,4 and 8. 
This higher oscillation amplitude will be seen to 
contrast with the lower A/D discussed above, in 
particular in the transition from an antisymmetric 
pattern of shedding to a symmetric. In line with 
our expectations, it is easier to observe symmet- 
ric shedding behind a rectangular cylinder rather 
than a square one. The symmetric modes obtained 
may be classified into three types, S-I to S-III. As 
mentioned earlier, the first two have been observed 
in experiments before, but on circular cylinders 
[HI [12] . The letter S signifies a symmetric pattern, 
while the number denotes how many pairs of shed 
vortices may be associated with one time period 
of the flow. Thus the cylinder sheds one vortex of 
each sign both at the top and the bottom in an S- 
II mode. At higher oscillation frequencies, the flow 
displays what we term as a mixed mode. 

To discuss the mixed mode, shown in figure [5j we 
choose a square geometry. Such a pattern has also 
been seen by Konstantinidis & Balabani [9 j. The 
shedding off the cylinder is actually symmetric, 
but some distance downstream, the shed vortices 
arrange themselves in an antisymmetric pattern, 
much like a Karman street, but with a larger spac- 
ing, and a correspondingly lower Strouhal number 
of 0.92 times that of a fixed square cylinder at this 




(d) fe/fo = 5 

FIG. 3: Vorticity fields at a typical time for A/D — 0.1 
at various excitation frequencies for a cylinder of aspect 
ratio 4. 
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FIG. 4: Flow patterns in the wake of an incline oscillat- 
ing rectangular cylinder at Re — 200 and A/D = 0.175. 
Circles: antisymmetric shedding, squares: symmetric 
shedding. The solid squares indicate the S-II mode, 
the open squares stand for the S-I mode, while the pat- 
terned square indicates an S-III shedding. Triangles: 
mixed mode, where the shedding is symmetric but the 
vortices arrange themselves into an antisymmetric pat- 
tern downstream. Stars: chaotic flow, single solid dia- 
mond: the Couder-Basdevant mode. 




FIG. 5: Mixed mode in the case of square cylinder, 
fe/fo = 4, A/D = 0.175. 



4, the wake pattern changes from a Karman street, 
followed by a chaotic pattern, through S-III and 
then S-II, followed by S-I with increase in the fre- 
quenc y of oscillation. The S-III mode, shown in 
figure 6(a)| is simply the S-II mode with an extra 
pair of vortices appearing close to the centreline. 
It is classified separately since it appears on the 
other side of chaos in the transition from antisym- 
metric shedding. A sample of S- II shedding on this 
cylinder is shown in figure |6(b)| 



The dominant frequency for f e /fo = 2 is f s = f Q 
(figurej7fc. This is indicative of subharmonic lock- 
on [6j [13] , while the shedding is locked on to the 
oscillation of the cylinder in the S-III mode. 



Reynolds number. The oscillating square cylinder 
together with the symmetric portion of its wake 
corresponds roughly to a stationary body with an 
effective D/h less than 1. In fact the Strouhal num- 
ber of the downstream portion of this figure is the 
same as that of a body whose aspect ratio is 0.67. 
In taller geometries, the mixed mode is actually en- 
couraged to occur by merger events of vortic es of 
one sign, some of which are evident in figures |3(b)| 



and |3(c) 



The downstream behaviour again be- 
comes antisymmetric. With all other parameters 
held constant, and reducing h alone, i.e., using a 
taller rectangular cylinder rather than a square, we 
would reduce the relative size of the boundary layer 
and therefore the strength of the shed vortex. The 
pressure oscillations, which normally promote anti- 
symmetric shedding, are correspondingly reduced, 
and so the symmetric pattern should persist fur- 
ther downstream for a given oscillation frequency. 
This is indeed manifested (not shown). 

Returning to our discussion on the aspect ratio of 



B. Mechanism for S-II mode 



Figure [8] shows the time signal of the vorticity u 
at the monitor point C. In the absence of a mean 
flow, it is easy to visualise the alternate shedding 
of oppositely signed vortices when the cylinder is 
moving to and fro. The mean flow advects both 
vortices downstream. This S-II mode of shedding 
is aided by the 'ground effect'. The primary vor- 
tices accelerate the fluid in the wake region towards 
the cylinder, and due to the larger area available 
on a rectangle rather than a square at the lee sur- 
face, significant vorticity of opposite sign is gener- 
ated. This effect is similar to the one studied by 
Carnevale et al [17] on a different problem. 

This ground effect is clearly visible on a rectan- 
gle of aspect ratio 8. The circles in black in fig- 
ure [8] indicate the time inst ances at wh ich vortic- 
ity field is plotted in figures 9(a)| - |9(d)| In figure 
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FIG. 6: Two of the modes of shedding at A/D = 0.175 
on a body of aspect ratio 4. (a) The S-III mode at 
fe/ fo — 2.15. Three pairs of binary vortices are shed, 
(b) The S-II mode at f e /fo = 4. In this mode two 
binary vortices are shed during each time period. 
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9(a) the cylinder is moving upstream. The pri- 



mary vortices are seen to form just behind the top 
and bottom surfaces of the cylinder. Vorticity is 
continuously supplied to them in the usual manner 
by the boundary layers. As these primary vortices 
grow they accelerate the fluid in the wake region 
leadi ng to the formation of boundary layers, figure 
|9(b)| on the lee side of the cylinder, of oppositely- 
signed vorticity with respect to the primary vor- 
tices. Now, when the cylinder moves downstream 
there is a local reverse flow near the cylinder. This 
causes shape changes in both vortices. The sec- 
ondary vortices continue to grow and cut off supply 
to t he pri mary vortices as can be seen in Figs. |9(c)| 
and |9(d)[ This cycle repeats. In the experiments of 
Xu et al [8] there was considerable overall reverse 
flow which aided the formation of opposite signed 
vortices on the surface of a circular cylinder. Here, 
aided by the ground effect, we obtain the S-II mode 
even without reverse flow at the inlet. 



FIG. 7: Power spectra at the monitor point A for 
D/h = 4 and A/D = 0.175. (a) Subharmonic shedding 
at fe/fo = 2. (b) The shedding is harmonic (symmet- 
ric) for f e /f = 2.15. 



C. Mode competition and chaos 



Perdikaris et al. PQ reported chaotic flow in the 
wake of a circular cylinder placed in a uniform flow 
at one particular amplitude of inline oscillation. In 
their simulations the cylinder was forced at the cor- 
responding Strouhal frequency of the fixed cylin- 
der. At about the same time, without being aware 
of that work, we had obtained chaotic flow using 
square and rectangular cylinders. We thus confirm 
their appealing finding. Further, while they sur- 
mised that a competition between antisymmetric 
and symmetric shedding was causing chaos, they 
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FIG. 8: Phase information for the S-II mode for D /h — 
8. Solid red line: the vorticity uj at monitor point B. 
Blue dashed line: inlet velocity U total- Here f e / f = 2 
and A/D = 0.175. The circles indicate the phases at 
which the vorticity field is shown in figures 9(a) - 9(d) 
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did not have direct evidence for this. In particular, 
their lift coefficients and spectra indicate antisym- 
metric shedding under all non-chaotic conditions, 
so they do not have mode competition between an- 
tisymmetric and symmetric modes. The use of a 
rectangular cross-section makes it easy to obtain 
symmetric shedding, so we are able to demonstrate 
that the shedding is antisymmetric at f e less than 
for the chaotic flow, and symmetric for f e greater 
than this value, which is a direct demonstration 
of mode competition in the sense of Ciliberto and 
Gollub [3H]. 



In figure 11(c), typical delay plots, for u y at 



fe/ fo = 2 and 2.085 are shown. The axes on the de- 
lay plots represent V\ = u y (t+r) and V2 = u y (t—r) 
at a s uitably chosen location and delay time r. In 
figure |10(a)| the paths in phase space are closed, 
which indicates periodicity. The noise in the com- 
putations gives rise to a patch rather than a single 
path, as often happens in these computations. In 
spite of this noise, this fi gure may easily be con- 
trasted with figure 10(b) , which is indicative of a 
chaotic flow. The chaotic window is easily visu- 
alised in figure [TT] with the antisymmetric mode, 
locked on to 0.5/ e , and the S-III mode, locked on to 
/ e , on either side of the narrow window of chaotic 
flow in between. 



Thus far we have obtained qualitative indications 



FIG. 9: The S-II mode of vortex shedding; D/h = 
8, f e /f = 2 & A/D = 0.175. Time has been 
non-dimensionalized using the convective time-scale, 
D/Uoo- (a) Attached primary vortices are growing, (b) 
Vorticity is generated on the lee side, (c) Close to the 
cylinder the flow is from right to left. The primary 
vortices are pushed apart and the secondary vortices 
are 'stretched', (d) The secondary vortices cut off the 
supply to primary. 



of the chaotic window. To confirm that the flow 
is indeed chaotic, we use a global measure S(At), 
defined as: 

S = ^^[uj(xi,yj,t + At) - uj(xi,yj,t )} 2 AxAy. 

i 3 

(2) 

Being an integrated quantity over the entire do- 
main, S is a reliable measure of chaos. For a peri- 
odic system of period T, S(T) = 0. Moreover for 
any At we should have S(At + T) = S(At). For 
fe/ fo — 2 and f e /f = 2.15, both these properties 
are seen in figure [l2j In particular, a sharp dip in 
S for At = nT for any integer n is visible. On the 
other hand, for f e /f Q = 2.085, S remains at a high 
value, characterising a chaotic system. 

A different behaviour is seen at D/h = 8. The 
vortex shedding mode changes from S-II to the 
mixed mode with S-I shedding when f e /f Q is var- 
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FIG. 10: Delay plots for D/h = 4, f e /fo = 2 and 
2.085. The oscillation amplitude A/D = 0.175. The 
monitor point is behind the cylinder at (0.5D, 0.35D). 
The delay plot (a) consists of closed curves, indicating 
periodicity, whereas (b) is characteristic of an aperiodic 
time signal. 



ied from 2 to 3.5. Figure 14 shows that the flow for 
fe/fo — 2 is periodic. However, the flow is chaotic 
for fe /fo = 2.5, as evidenced by the spectrum in 
figure |13(a)| Also the arrangement of vortices at 
a given time is in no particular pattern, as seen 
in figure |T3j The contrast in terms of the S is 



demonstrated in figure 14 A smaller region is cho- 
sen to improve the contrast. The variation of S 
for fe/fo = 2 is seen to be periodic, whereas for 
fe/fo = 2.5 it is not. 



D. Conclusions 



To summarise, we have studied two-dimensional 
flow past inline oscillating cylinders of rectangular 
cross-section at a relatively low Reynolds number. 
The S-II mode of shedding is obtained computa- 
tionally for the first time to our knowledge. The 
S-II mode is enhanced in the case of rectangular 
cylinders due to the ground effect. A variant sym- 
metric mode, named here as S-III is also observed. 
The type of shedding changes as the frequency of 
cylinder oscillation is increased, all other parame- 
ters held constant. At lower oscillation amplitudes 
the frequency regime in between the antisymmetric 
and symmetric modes consists of a periodic flow 
where the shedding is neither symmetric nor an- 
tisymmetric, but a constant phase is maintained 
between shedding on the upper and lower surfaces. 
At higher oscillation frequencies, the flow in this in- 
termediate frequency regime is chaotic. More than 
one window of chaos may exist in the frequency 
range. Since these windows always lie between 




(b) 




(c) 

FIG. 11: Chaotic window for D/h - 
the shedding here is antisymmetric, 
the shedding is chaotic, (c) f e /fo : 
of symmetric shedding. 



= 4. (a) f e /f = 2, 
(b) fe/fo = 2.085, 
= 2.15, S-III mode 



regimes of antisymmetric and symmetric shedding, 
the chaos is due to mode competition in the sense 
of Ciliberto & Gollub [3 j. The shedding frequency 
in the periodic regimes on either side of the chaotic 
window are locked-on to two different submultiples 
of the excitation frequency. 



The use of a rectangular cylinder rather than a 
square one has made several of the above obser- 
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vations possible. A global, and therefore reliable, 
measure has been used to ascertain the existence 
of chaos. It is hoped that the present results will 
motivate experiments with rectangular oscillating 
cylinders. The work also presents natural exten- 
sions to flow past accelerating bodies. 



FIG. 12: D/h = 4. Variation of S with At for different 
cases. The black curve is for f e /fo = 2, the red curve 
for fe/fo — 2.085 and the blue curve for f e /fo = 2.15. 
The black and blue curves exhibit periodicity and have 
sharp dips at At = nT for any integer n, whereas the 
red curve is aperiodic. 
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FIG. 13: D/h = 8, f e /f = 2.5 & A/D = 0.175. 
(a) The spectrum is broadband, and (b) the arrange- 
ment of vortices is unordered, showing that the flow is 
chaotic. 



FIG. 14: D/h = 8. Variation of S with At in a region 
x/D = 2.5 - 8.75. The blue curve is for f e /f Q = 2, and 
is repetitive indicating periodicity. The red curve is for 
the case f e /fo = 2.5, and shows the aperiodic nature 
of the vorticity field in that region. 
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